NASA TECHNICAL NOTE 


NASA TN 0-7666 



N 74- 2 1290 

Unclas 

36622 


LANDING RATES FOR MIXED 
STOL AND CTOL TRAFFIC 


ffC&h S A -T N - D - 7 6 6 6 ) LANDING HATES'" FOR MIXED 

STOL AND CTOL TBAFFIC (NASA) 76 p HC 
$4-00 - CSCL 17G ~ 

. H 1/2 1 


by John S. White 

Ames Research Center 
Moffett Field, Calif. 94035 







NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • 


WASHINGTON, D. C. 


APRIL 1974 


1. Report No. 


2. Government Accession No. 


3. Recipient's Catalog No. 


D- 7666 



4. Title and Subtitle 

5. Report Date 


APRIL 197 h 

LANDING RATES FOR MIXED STOL 

AND CTOL TRAFFIC 

6. Performing Organization Code 

7. Author(s) 

8. Performing Organization Report No. 

John S. White 


A-5044 


10. Work Unit No. 

9. Performing Organization Name and Address 


NASA Ames Research Center 

501-03-11 

Moffett Field, Calif. 94035 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Technical Note 

National Aeronautics and Space Administration 

14. Sponsoring Agency Code 

Washington, D. C. 20546 



15. Supplementary Notes 


16. Abstract 

A study was made to determine the expected landing rate for STOL-only 
traffic and mixed STOL-CTOL traffic. The conditions used vary from present 
day standards to an optimistic estimate of possible 1985 conditions. A 
computer program was used to determine the maximum landing rate for the 
specified conditions and aircraft mix. The results show that the addition 
of STOL on a CTOl runway increases the total landing rate if the STOL airborne 
spacing can be reduced by use of improved navigation equipment. Further, if 
both takeoff and landings are performed on the same runway, the addition of 
STOL traffic will allow an increase in the total operation rate, even with 
existing spacing requirements. 


17. Key Words (Suggested by Author(s)) 

Air Traffic Control 
Terminal Area Guidance 
Aircraft Operations 


18. Distribution Statement 

Unclassified — Unlimited 

CAT. 21 

19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 


76 

$4.00 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 

























TABLE OF CONTENTS 


Page 

NOTATION v 

SUMMARY 1 

INTRODUCTION 1 

THEORY OF LANDING-RATE CALCULATIONS 3 

Airborne Spacing 3 

Runway Clearance 6 

Probability of Abort 8 

Landing Rate 8 

Total Operation Rate 9 

COMPUTER PROGRAM : 9 

PROGRAM INPUTS 10 

RESULTS 12 

Single-Velocity Approach 14 

Dual-Velocity Approach 17 

CONCLUSIONS 20 

APPENDIX A - DERIVATION OF EQUATIONS FOR T l2 22 

APPENDIX B - PROGRAM LISTING AND DESCRIPTION 48 

REFERENCES 70 


PRECEDING PAGE BLANK NOT FILMED 


A-5044 iii 



NOTATION 


Fortran Name Symbol 


ACC 

a 

DA 

Da 

SIG 

K 

M or ZM 

m 


N lz 

RD . . . 

r(-) 


S 

SF 

S f 

SO 

So 


t 


T 


T a 

T2ACT 

T act 

To 

To 


T 

m 

a 


T 

m 

TNOM 

T 

nom 

TTO 

T, 

to 

TVF 

T , 

vf 

T12A 

T\2 

T12M (I , J) 

Til 

nom 

VA or VB 

V 


aircraft acceleration capability, m/sec 2 

distance from runway threshold at which velocity is 
changed, km 

constant used to calculate AT 

distance from runway threshold to outer marker, km 
number of takeoffs between successive landings 
range or expected variation of stochastic quantity 
required separation distance, km 
required separation at touchdown, km 
required separation at outer marker, km 
time, sec 

average time between landings, sec 

time before touchdown at which velocity is changed, 
sec 

actual time from outer marker to touchdown, sec 
runway clearance time, sec 
time when x\ = m 2 , sec 
time when x = m, sec 

nominal time from outer marker to touchdown, sec 
time required for takeoff, sec 

time at which aircraft must be at final velocity, sec 
actual interarrival time, sec 
nominal interarrival time, sec 
aircraft velocity, m/sec 


PRECEDING PAGE BLANK NOT FILMED 


v 



X 

aircraft location, km ' 


X ' 

location of separation boundary, km 

DELP 

6p 

error in position at time of nominal passage over 
outer marker, km 

DELT 1) 
DELT 2) 

6T 

landing time deviation, sec 

DEL. . 

«(’) 

error in ( ) 

TDL 1 1 
TDL 2 J 

AT 

allowance for errors in landing time, sec 

ZLR 

■ h 

landing rate, aircraft/hr 

ZNOPR ' 

X 

o 

operation rate (landings plus takeoffs), aircraft/hr 


e 

generalized error 

TDEL 

T 

C 

interarrival time required for runway clearance, sec 

TFIN 

1 fin 

interarrival time required for airborne spacing at 
touchdown, sec 

TIN 

T . 

,%n 

interarrival time required for airborne spacing at 
merge point, sec 

TMED 

. T med 

interarrival time required for airborne spacing at 
velocity change point, sec 

Subscripts 



1,2 


refers to first or second aircraft in the landing 


sequence 


vi 


LANDING RATES FOR MIXED STOL AND .CTOL TRAFFIC 


John S. White 
Ames Research Center 


SUMMARY 

A study was made to determine the expected landing rate for STOL-only 
traffic and mixed STOL-CTOL traffic. The conditions used vary from present 
day standards to an optimistic estimate of possible 1985 conditions. A com- 
puter program was used to determine the maximum landing rate for the specified 
conditions and aircraft mix. The results show that the addition of STOL on 
a CTOL runway increases the total landing rate if the STOL airborne spacing 
can be reduced by use of. improved navigation equipment. Further, if both 
takeoff and landings are performed on the same runway, the addition of STOL 
traffic will allow an increase in the total operation rate, even with existing 
spacing requirements. 


INTRODUCTION 


Numerous studies have been made to determine the landing rate of aircraft 
on a single runway. In one of the early studies, Blumstein (ref. 1) con- 
sidered a traffic mix of various commercial transports arriving at an airport 
and computed the landing rate, assuming that aircraft were always available at 
the gate. Simpson (ref. 2) and Odoni (ref. 3) describe theoretical equations 
for landing rate, including the effects of terminal area operations beyond the 
gate, and assume that the aircraft arrivals have a Poisson distribution. Both 
give some numerical results, but for CTOL vehicles only. An additional study 
by the National Bureau of Standards (ref. 4) used a different concept of land- 
ing capacity, which involved a facility serving sequential customers of vari- 
ous types. Again, however, the results presented are for CTOL vehicles only. 
The capacity of a STOL runway was considered by Rinker (ref. 5), but numerical 
results were presented for only two cases. 

The purpose of this report is to supply additional information for STOL- 
only runways and to consider the effect of mixed STOL-CTOL traffic. Although 
the National Aviation System Plan (ref. 6) implies that STOL and CTOL traffic 
will use separate runways, there may be times when a common runway is required 
on a temporary basis. Thus this report develops a general theory of landing 
rates, which can be applied to separated or mixed traffic, and from which 
expected landing rates can be determined. Further, it indicates the accura- 
cies in control and navigation required to achieve various levels of landing 
rates. Finally, previous studies have assumed that all aircraft fly at con- 
stant airspeed all the way from the instrument landing system (ILS) gate to 
touchdown, a distance of 6 to 10 miles. An alternate approach is for an air- 
craft to fly at 160 to 180 knots at the ILS gate and then reduce speed gradu- 
ally to the final touchdown speed. The effect of such a speed change on the 
landing rate is also considered. 



The general approach used here parallels that of Blumstein, but it is 
considerably more complex to account for the variable speed approach and dif- 
ferent outer markers for CTOL and STOL. Also, Blumstein had only two separa- 
tion criteria — airborne spacing at the outer marker and runway clearance. 

This study adds another criterion: adequate airborne spacing must be main- 

tained until touchdown. 

In this study, it is assumed that landings have priority, and that the 
aircraft land on a "first come, first serve" basis. We are concerned with 
maximum capacity, and assume that the aircraft are immediately available to 
start a landing approach whenever desired. This approximates the situation 
in which a good 4-D guidance system controls the aircraft in the terminal area 
to arrive within a few seconds of the desired time. If these arrival errors 
become large, they will approach a Poisson distribution, and the landing 
capacity will be reduced. 

Takeoffs are also considered, but only when the sequence of landings is 
such that a takeoff can be inserted without affecting the landing rate. 

Sketch (a) shows some possible approach paths. The various CTOL 
approaches merge fairly far out at the CTOL merge point, while STOL vehicles 
make curved approaches and arrive on the runway centerline at the STOL merge 
point. 

The results of this study indicate the effect on airport capacity of 
varying the separation requirements and navigation accuracy. To fully utilize 

the increased runway capacity sug- 
gested herein will probably require a 
great deal of additional automation, 
both on the ground and in the air. 

This automation will probably not be 
available before about 1985. Further, 
in determining the actual separation 
criteria, additional factors, such as 
trailing vortices, pilot and control- 
ler reaction times, accuracy of Air 
Traffic Control (ATC) surveillance 
radars must be considered. However, 
the optimistic viewpoint taken here 
does provide useful information to 
aid in the determination of desirable 
trends in separation requirements, and 
indicates what portion of the approach 
system is presently acting as the 
limiting factor. 

The remainder of this report 
first develops the theory of landing- 
rate calculations, then briefly 

Runway 


\ 

\ 



Sketch (a) 
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describes a computer program to solve these equations, and presents the 
inputs and results obtained. 


THEORY OF LANDING-RATE CALCULATIONS 


To determine landing rates, some of the characteristics of a single air- 
craft approaching the runway are first considered, then a technique is devel- 
oped for determining the minimum spacing between successive aircraft, and, 
finally, the maximum landing rate associated with a specific mix of aircraft 
is determined. 


Airborne Spacing 



may occur anywhere along this velocity time profile; two possible locations 
are shown in the figure. 

To simplify the equations associated with changing velocities, it was 
assumed that the aircraft would fly at constant speed V a until it was a dis- 
tance D a from the end of the runway, at which time it would instantaneously 
change velocity to l/fr. This step change in velocity simplifies the calcula- 
tion and, for reasonable values of separation, decreases the separation by 
less than 10 percent. 

To determine D a , first locate the time T a at which the two constant- 
velocity segments intersect, that is, T a = T v f + ( V a - V^lla. Since 

D a = T a Vb (fig. 1), 


D a = T v f v a 



CD 
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The location of the aircraft, x, at any time t is specified by 



f-v 

> 

x ± D a j 

X = • 

I- (■ ■ 1 ) 

~ V aP > 

x > D a J 


( 2 ) 


T = 
m 


time 

the 

aircraft is 

over 

the 

m 

' V b 


> 

m 

< D a 

m 


(i - M , 

m 

> D„ 

V + 
v a 

v a 1 

^ V 

a 


( 3 ) 


Consider two aircraft approaching a given runway in sequence. A minimum 
separation distance is required between the two aircraft at all times during 
the approach and, further, the second aircraft must not touch down until the 
first aircraft is off the runway. The problem is to determine the minimum 
time between touchdowns that meets the separation criteria and, from this, to 
determine the maximum landing rate. The two aircraft may be of different 
types, with different approach paths, velocities, time-to-clear runway, sep- 
aration distances, etc. The separation requirements are stated in terms of 
in-trail distances and must be met only while the two aircraft are on a 
common path; before that, the vehicles are considered to have satisfactory 
lateral separation. - 


A few terms should be defined more explicitly.. . First, the merge point 
is that point where the aircraft arrives on the extended; runway centerline, 
and is. stabilized on the runway heading. The . location-, -of the merge point, m,.. 
varies. from one type of aircraft to another. The common path is the section, 
of the flight path that will be used by two successive aircraft; therefore, .. 
it occurs between the minimum of and m 2 and the runway and its length 
depends on the types of the two aircraft. ...... 


The separation requirements are defined so that, while- the first aircraft 
is on the common path, the second. aircraft must be far enough behind. The 
present-day separation requirement during approach is 6j, km (~3 miles). Since 
the navigation systems of STOL vehicles may be more accurate, we wish to con- 
sider the possibility of reduced distance separations at the merge point, and 
also a smaller separation requirement at touchdown than at the merge point. 

Separation is defined as that portion of airspace reserved for a given 
aircraft and is assumed to be half before and half behind the nominal air- 
craft position. The separation criterion for two aircraft is shown in 
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figure 2, and the requirement. is 
that the reserved airspaces do not 
overlap, or 

= xi + S\/2 < x '2 = x 2 - S 2 /2 

Figure 2.— Separation criteria. 

where S\ and S 2 , the separation 
distances associated with each 

aircraft, need not be constant with time. For this study, define a separa- 
tion requirement, So^> at the merge point, and Sfo, at touchdown, and assume 
that it varies linearly between those two points, that is. 


Reserved oirspoce 
A/C I 


Reserved oirspoce 
A/C 2 


S,/2 


(4) 




V 



for x . < m . 
^ ^ 


x . > m . 


(5) 


This form for £ = 1,2 was chosen to 
allow for a possible reduction in 
spacing as the aircraft approaches 
the threshold while maintaining rel- 
atively simple equations. 

To visualize the separation 
requirement, figure 3(a) shows a 
specific example of two approaching 
aircraft. For each aircraft, its 
merge point and velocity change 
point are shown. To guarantee sep- 
aration at all times, checks must be 
made at three places, namely, at the 
end of the common path, at the begin- 
ning, and at the point where the . 
velocity of aircraft 2 changes. In 
the example shown, these tests occur 
at t = 0, t = 2V (since the common 
path starts at tne closest merge 
point, m 1 , as drawn), and at t = T a ^. 
At these points, we want to guarantee 



(a) Generalized situation. 

Figure 3.— Separation between 
that successive aircraft. 


a: 2 - = r 2 - 


- ^1 



( 6 ) 


For different types of approaching aircraft, the relative locations of 
m-i , mo, D , and D , vary, as does the time of the tests, resulting in a 

1 ^ <Z\ o. 2 

variety of situations, each of which must be considered individually. 

A-5044 
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To determine the closest possible spacing, three interarrival times are 
calculated by selecting x 2 , in inequality (6), so that the difference is zero 
at each of the three test times. The largest of these times guarantees ade- 
quate separation throughout the airborne portion of the flight. The effect 

of runway clearance requirements is 
considered in the next section. 




Considering the example of 
figure 3(a) again, first calculate 
t fin, which assumes that minimum sep- 
aration occurs at t = 0. This is 
shown in figure 3(b), where the tra- 
jectory for aircraft 2 is moved 
closer to that for aircraft 1, and 
the minimum separation, {Si + S 2 )l2, 
occurs at t = 0. Next, calculate 
t i n , which assumes minimum separa- 
tion at t = T mi (fig. 3(c)). Finally, 
calculate x me ^, which assumes mini- 
mum separation at t = T ai (fig. 3(d)). 
For this example, the critical test 
is the calculation for x„. , and an 
interarrival time of T 12 - T fi n > 
guarantees adequate airborne^ separa- 
tion throughout the flight. 


The various cases that can 
occur, depending on the locations of 
m T>~m jr r ~Dn i > and D n 2 T ~ar e~ c le scribed 


in appendix A, along with the equa- 
tions for calculating the various) t. 


(c) Minimum separation at 




(d) 


Minimum separation at 



Figure 3.— Concluded. 


Runway Clearance 

Having computed the t value that 
guarantees satisfactory separation 
everywhere on the flight path, we 
now consider separation on the 
runway. The FAA rule is that a fol- 
lowing aircraft may not touch down 
until the preceding aircraft clears 
the runway. The time a given type 
of aircraft needs on the runway 
before turning off, T c , is known, 
but, in scheduling aircraft arrivals, 
some time must be added to allow for 
the possibility of a late landing of 
the first aircraft or an early land- 
ing of the second aircraft. Four 
sources of error in landing time are 
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considered: (1) a position error at the merge point, 6p, (2) an error in the 

approach velocity, 6V a , (3) an error in the touchdown velocity, 67£, and (4) 
an error in the time at which the aircraft reaches the final velocity, ST v f. 
For a given aircraft, the nominal time, T nom , from its merge point to touch- 
down is given by the negative of equation (3) , and the actual time is 


^aot < 


m + Sp 


m + 5p D a + & D a 


m < D a 


V a * 6V a V a + 6F a 


1 - 


^a + 

V, + 67, 
b b 


m > D. 


a 


(7) 


where 


Da + ™a = W vf + *T vf ) Wj, + « V + 


[(V a + 67 a ) - (V b + 67 fc )](7 fc + 67 fc ) 


2a 


which results in a landing time error of 


^ ~ ^act ~ 


'-nom 


( 8 ) 


From this, the actual time between touchdowns can be calculated, given the 
nominal time, by 


12 


act 


= T 


12 


nom 


- 6T l + 


6T. 


(9) 


and if ^l 2 a at •*‘ ess t ^an the runway clearance time, T c , the second aircraft 
must abort its landing. This implies that an allowance must be made for <5 Tj 
and ST 2 to minimize aborts. Assuming the four error sources are independent, 
one can take the partials of T\2 aQ t with respect to the various error sources, 
multiply them by the expected range of the error, and add them on an RSS basis 
to obtain a time allowance, A T: 




(10) 


where 0 represents one of the four error sources, rg^ is the range of prob- 
ability distribution of that error source, and K is chosen so that A T will be 
large enough to minimize the probability of abort. The planned separation 
for runway clearance, t , should then be 
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t c = T e + Ar (11) 

and the nominal interarrival time can be obtained from the various x values 
as 


T\ 2 nom max ( x in, x med, x fin, x c ) (12) 

where Ti 2 M depends on the types of aircraft 1 and 2. (The equations for 
the various partials used in equations (10) are given in appendix A.) 


Probability of Abort 

Since we are trying to maximize the number of aircraft to be landed, the 
aircraft will be flying close together, and if one separation criterion is 
violated, the second aircraft should abort the landing. It is assumed that 
airborne separation has been determined from a consideration of the likely 
errors in position and velocity of the aircraft; therefore, aborts due to 
violation of the airborne separation will not occur (at least from the 
expected errors) . However, the runway occupancy rule does not consider pos- 
sible errors; any violation of this rule will therefore cause an abort, as 
implied by equation (9). The quantity AT (eq. (11)) is designed to consider 
the effect of errors on the runway occupancy rule, and K (eq. (10)) is chosen 
so that the probability of abort because of a violation of the runway occu- 
pancy rule will be sufficiently small. As shown in appendix A, if K- 1, the 
probability of abort is less than 0.006. 


Landing Rate 

To determine the average interarrival time T, calculate 1\z norn for all 
possible sequences of aircraft types, and then obtain an average, weighted 
according to the relative probabilities of arrival associated with each type 
of aircraft. It is assumed that there are n types of aircraft to be con- 
sidered, and the probability that the next aircraft is of type i is P^ and is 
independent of the type of the previous aircraft. One then obtains 


T 


n 


n 


EE 

i=l J=l 


P.P.T.. 

1 <7 T'C 


nom 


(13) 


where obtained from equation (12), depends on the types of aircraft 

i and j, and n 



i=l 


The nominal maximum landing rate, A , is then 


\ L = 3600 /T 


8 


(14) 



which assumes that aircraft are continually available to start their approach 
as soon as needed, and is thus a true measure of runway capacity. To achieve 
this capacity the terminal area controller must- be told when an aircraft-, 
should be at the outer marker, and a time-controlled guidance scheme, such as 
suggested by Erzberger and Pecsvaradi (ref. 7), must be available. 


Total Operation Rate 

Another factor of interest is the total operation rate, X 0 , of a runway 
which allows takeoffs and landings to be mixed. To determine X 0 , consider a 
specific interarrival time and determine if it is long enough to allow a 
takeoff to occur without delaying the landing of the second aircraft (there- 
fore, landings have priority over takeoffs). Assume that the time for take- 
off is constant? for all aircraft, and that the types of aircraft 1. and 2 
are specified. We can determine the largest integer that satisfies 

-■ Nl2 T t o < C^l 2- nom ~ T a ) (15) 

which specifies that aircraft can take off in the available interarrival 
time. 1 From the weighted sum of over all types of aircraft, as in 
equation (13), and the landing rate, X^, the total operation rate, A 0 , is 
obtained: ' 

■ X L f 1 * E F i P /ij) C16) 

t- , j 

COMPUTER PROGRAM 


A computer program was written to solve the various equations given 
above so that the landing rate and total operation rate could be determined 
as a function of the various parameters of the system. The program is writ- 
ten in Fortran IV, and a brief description and listing are included in 
appendix B. 

For each type of aircraft (the program handles up to 12 types), the 
following input quantities are specified: 

m distance of the outer marker from runway threshold 

So,Sf airspace required for separation at the outer marker and at 

touchdown 


^his takeoff criterion is not strictly compatible with the present FAA 
requirements that specify 1 minute between successive takeoffs and 2 miles 
separation between a takeoff and the following landing. However, it suggests 
how arrivals can be interspersed with landings. 
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T. 


runway clearance time 

approach velocity and touchdown velocity 

time before touchdown when the aircraft should be at velocity V ^ 


Q 

Va,V b 

V 

a Irange of uncertainties of the various quantities 
P probability of arrival of that aircraft type 


The principal outputs were calculated from deterministic data using weighted 
averaging techniques. These outputs were Xjj, the landing rate, and \ a , the 
operations rate (calculated for T-^ 0 = 36 and 48 sec) . Additional output was 
obtained from Monte Carlo type runs in which 300 landings were made to deter- 
mine the number of landing aborts that occurred for a given set of conditions 
For these runs, the random variables Sp, &V a , etc., were selected from a tri- 
angular distribution with the width equal to the corresponding range. Also, 
to control the number of aborts, the value K in equation (10) was variable. 


PROGRAM INPUTS 


To obtain average landing rates, the program referred to in the previous 
section was used over a range of the various input parameters. The specific 
values used, along with a justification for the choice, are indicated in this 
section. The parameters are the aircraft mix, and the quantities m, S 0 , Sf, 
T c , v a> v b> T vf> a > r *>P> *&Va, v&Vb* and r $tvf selected separately for each 
type of aircraft. 

For CTOL vehicles, two values of distance to the merge point, m, were 
used — 19 km (=10 miles) and 7 km (=4 miles), which represent fairly long 
and fairly short approach paths presently in use. For STOL vehicles, 7 km 
(=4 miles) and 2 km (=1 mile) were used since it is unlikely that STOL vehi- 
cles will use the long approaches of the CTOL vehicles and, when a curved 
approach is used may turn onto the final heading fairly close to the runway. 

The nominal values of separation, S Q , are 6 km (*3 miles) for CTOL vehi- 
cles and 2 km (*1 mile) for STOL vehicles. The CTOL separation is represen- 
tative of present standards, but STOL vehicles may have more precise naviga- 
tion and control systems so that reduced separations might be possible. 

The airborne separation at touchdown, Sf, was usually set equal to S 0 , 
the separation at the outer marker, but some runs were made with Sf essenti- 
ally zero. In this case, the actual separation at touchdown is held to rea- 
sonable values by both S Q and T a . It is not clear whether present practice 
is to have Sf = S 0 , or to ignore it. 

The runway clearance time was determined by computing braking time and 
adding to that an allowance for taxi time to the turn off. Braking time to 
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a dead stop, using 0.25 g deceleration, is roughly 30 sec for CTOL and 15 sec 
for STOL. The assumption of an additional 30 and 15 sec, respectively, for 
taxi time gives the runway clearance time used, namely, 60 sec for CTOL vehi- 
cles and 30 sec for STOL. Data were also obtained for high-speed turnoffs 
that are immediately available giving runway clearance times of 30 sec for 
CTOL and 15 sec for STOL. 


Two different situations were assumed for the approach velocity. First, 
the aircraft is assumed to fly at touchdown velocity from its merge point 
and, second, the aircraft is assumed to be at a velocity commanded by ATC 
when it passes its merge point, and to slow down during the approach. The 
touchdown speeds, V ft, used were 60, 65, and 70 m/sec (=120, 130, and 140 knots) 
for CTOL and 30, 35, and 40 m/sec (-60, 70, and 80 knots) for STOL. It is 
felt that these represent present-day CTOL touchdown speeds and proposed STOL 
touchdown speeds. Also, it was assumed that, for the dual-velocity case, the 
initial approach velocity, V a , was 90 m/sec (sl80 knots) for all CTOL vehicles 
and 60 m/sec (=120 knots) for all STOL vehicles, which represent terminal 
area cruise speeds. 


In any case, a CTOL vehicle was required to reach its touchdown velocity 
1-1/2 min before touchdown and a STOL vehicle, 1 min before touchdown. These 
are felt to be the current requirements set forth by pilots, and were the 
values used for T v f. The deceleration, a , used during the change of speed is 
0.6 m/sec 2 (=>0.06 g) for all vehicles. 


As mentioned earlier, all random variables are assumed to be independent 
and to have a triangular distribution, centered about the nominal value, with 
a specified range. Thus, for a generalized random variable, 0: 


e 

nom 


re / 2 < e 


< 

act ~ 


0 

nom 


+ re/2 


It is assumed that the pilot will maintain velocity within ±1-1/4 m/sec 
(±2.5 knots) so that the range for both 5V a and is r6V a = r&V^ = 2.5 m/sec. 
The position at the merge point is assumed to be ±1 km (±1/2 mile), rSp = 2 km, 
and the variation on T v £ is assumed to be r&T v f = 30 sec for CTOL vehicles 
and 10 sec for STOL vehicles. 


A brief study of the actual mix of CTOL traffic at the San Francisco 
Airport was the basis for the CTOL traffic mix used here. Three different 
landing speeds were assumed and their relative percentages are shown in 
table 1. The same percentages were used for STOL traffic. For mixed STOL 
and CTOL traffic, the percentage of aircraft with each landing speed is shown 
in table 1. The mix of aircraft shown is fairly representative and all 
results presented here use this mix. Data were obtained for a mix with some- 
what higher average landing speed, which resulted in a slight increase in 
landing rate, but the trends were not affected. 
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TABLE 1.- AIRCRAFT MIX FOR VARIOUS PERCENTAGES OF STOL TRAFFIC 


Type 

STOL, 

^\oercent 
Landing^"-*, 
speed, m/sec^\ 

0 

25 

' 

50 

75 

100 

) 

30 

0 

11 

22 

34 

45 

STOL } 

35 

0 

9 

17 

25 

34 

J 

40 

0 

5 

11 

16 

21 


60 

45 

34 

22 

11 

0 

CTOL 

65 

34 

25 

17 

9 

0 


70 

21 

16 

11 

5 

0 


RESULTS 


The computer program was used with the input data just described for a 
number of cases. The results are presented in a series of plots, following a 
brief discussion of the effects of some of the uncertainties. 

The uncertainties in position and velocity affect the landing rate only 
through the calculations to guarantee adequate spacing on the runway (air- 
borne spacing was controlled by the required separation). As a result, 
reducing these uncertainties had an effect very similar to reducing the time 
on the runway. 

The nominal position uncertainty of 2 km at the merge point increased 
the required time spacing between touchdowns by between 15 and 30 sec, 
whereas each of the other uncertainties increased it by about 1 or 2 sec. If 
these numbers are compared with the nominal runway occupancy time of 30 to 
60 sec, it is apparent that reducing the position uncertainty (or, equiva- 
lently, increasing the accuracy of time of arrival at the merge point) might 
change the landing rate significantly, but reducing the other uncertainties 
would have virtually no effect. Therefore, only the results for position 
uncertainty are presented here. 

The basic results are presented in figures 4 through 9 for single-velocity 
cases and in figures 10 through 15 for dual-velocity cases. In each figure, 
the landing rate is plotted against the percentage of STOL traffic. For sim- 
plicity, only two landing rate curves are shown which represent merge point 
distances of 19 km for CTOL and 7 km for STOL, and 7 km and 2 km, respec- 
tively. Other combinations of merge point distances fall between the curves 
given. 

The figures also show the total operation rate (sum of takeoffs and 
landings) for takeoff times, Tf, 0> of 36 and 48 sec. These times approximate 
minimum and maximum takeoffs for CTOL jets- and, although a STOL vehicle could 
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take off in less time, it was felt that an allowance of less than 36 sec for 
takeoff would be unrealistic. These curves are shown as bands in the figures 
rather than as individual curves 
because of the step-wise fashion 
in which takeoffs are added. This 
procedure is explained in sketch (b) 
where landing and total operation 
rate are plotted versus separation 
for two values of m. The conditions 
for sketch (b) are: two types of 

STOL vehicles, each using a single- 
velocity approach, one at 40 m/sec 
and the other at 30 m/sec, with 
50 percent probability of each type, 

T c = 30, Tfc 0 =36. As the landing 
rate decreases, because of an 
increase in separation, the number 
of takeoffs that can be inserted 
between landings increases in a 
step-wise fashion, one takeoff at a time. These steps occur at different 
values of S Q for m = 2 and m = 7 km, so that the total operation rate is 
sometimes larger for one case and sometimes for the other. Other values of 
m would cause the steps to change at other locations. As a result, it is 
very difficult to determine the effect of m on the total operation rate and, 
as indicated previously, the total operation is shown as a band in figures 4 
through 15. One noticeable trend is that the total operation rate tends to 
increase as the landing rate decreases. 

Table 2 presents the values of the parameters used for each case. Six 
cases are listed and each was run for both the single-velocity and dual- 
velocity approach. 


TABLE 2.- DATA FOR EACH CASE 


Case 

STOL 
S 0 , km 

CTOL 
Sq y km 

STOL 
T a , sec 

CTOL 
T c , sec 

STOL 
r&p, km 

CTOL 
r6p, km 

Present 

6 

6 

30 

60 

2 

2 

Nominal 

2 

6 

30 

60 

2 

2 

Reduce runway 
time 

2 

6 

15 

30 

2 

2 

Improve STOL 
navigation 

2 

6 

30 

60 

1 

2 

Minimal navi- 
gation errors 

2 

6 

30 

60 

.1 

.1 

Minimal spacing 

.1 

.1 

30 

60 

.1 

.1 



S 0 „ Km 

Sketch (b) 


The present case used the current FAA rules for both STOL and CTOL. For 
the nominal case, since it is assumed that the STOL vehicles will use more 
accurate path guidance, the STOL spacing can be reduced to 2 km. Next is a 
nominal case with the runway occupancy time cut in half, then the case for. 
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which the STOL navigation and guidance errors were reduced from the nominal. 

To show the maximum effect of these errors, they were essentially eliminated 
for both CTOL and STOL for the next case. For the final case, with the spac- 
ing and errors reduced to nearly zero, the aircraft are not allowed to pass 
and must obey the runway occupancy rule, but otherwise they may be extremely 
close. Separation requirements for trailing vortices are not considered. 

The results are now discussed, case by case, to show the changes to the 
landing rate and total operation rate, as STOL traffic is added to existing 
CTOL traffic. Also, the factor that controls the landing rate will be 

described, that is, whether the con- 
trolling factor is the airborne spac- 
ing at the merge point, at the touch- 
down point, or, for the dual-velocity 
case, at the midpoint (where the 
velocity changes), or whether it is 
the runway occupancy time. 

Single-Velocity Approach 

Present standards— The situation 
using present standards is shown in 
figure 4. As expected, increasing 
the percentage of STOL traffic reduces 
the landing rate even when the merge 
point for CTOL and STOL traffic is 
fairly close in. Also, for fixed 
percentages of STOL traffic, reducing 
m STOL an d ^croL increases the landing 
rate. The factors that control the 
landing rate are airborne spacing at 
the merge point when the second air- 
craft is slower than the first and 
the spacing at touchdown when the 
second aircraft is faster. In either 
case, the spacing is such that runway 
occupancy time does not affect the 
landing rate. 

When total operations are considered, the addition of STOL traffic may 
allow an increase in the rate, assuming aircraft are waiting to take off, 
since there are more gaps in the arriving traffic which can be used for 
takeoffs. 

Nominal— The nominal case, with improved STOL navigation accuracy and 
therefore reduced STOL separation, is shown in figure 5. Here, adding STOL 
traffic decreases the landing rate only slightly and, with the higher percen- 
tages of STOL traffic, the landing rate actually increases. Also, an all- 
STOL runway has a higher landing rate than a CTOL runway. The controlling 
factor is again the spacing at the merge point or at touchdown when the first 
aircraft is STOL. 



STOL, percent 

Figure 4. — Single velocity, 
present standards. 
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When total operations are considered, the addition of STOL vehicles 
increases the total operation rate. 

Reduced runway time— Next consider the case with reduced time on the 
runway for both CTOL and STOL (fig. 6). Since the landing rate for CTOL vehi- 
cles is controlled by airborne spacing, reducing the runway time does not 
affect the landing rate although it allows a considerable increase in the 
total operations. The landing rate for STOL vehicles was controlled by runway 
time so that the STOL landing rate increased considerably. A further reduc- 
tion in the runway time below 15 sec would increase landing rate only slightly. 
Again, for this case, introducing STOL vehicles to the traffic mix decreases 
the landing rate only slightly, with a considerable increase for larger per- 
centages of STOL. The trend is for increasing total operations as STOL air- 
craft are added to the mix. 



Figure 5. — Single velocity 
nominal case. 


Figure 6.— Single velocity reduced 
runway time. 
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Improved STOL navigation— In figure 7, the nominal case is considered 
with improved guidance and navigation for the STOL vehicles, which reduces 
the initial position error from 2 to 1 km, with a corresponding reduction in 
the uncertainty of the landing time. Since this uncertainty is directly added 
to runway time (eq. 11), the landing rates here and for the previous case 
(fig. 6) are essentially the same. The effect on total operations is quite 
different, however. For CTOL only, the total operation rate is the same as 
for the nominal case (fig. 5), and the increase with increasing STOL traffic 
is similar to that for the nominal case. 

Minimal navigation error— In figure 8, the navigation errors for both 
STOL and CTOL aircraft were set to essentially zero. The resulting landing 
rate and total operation rate are essentially the same as for the previous 
case (fig. 7), indicating that an extreme reduction in the navigation error 
alone is not useful. 



Figure 7.— Single velocity improved Figure 8.~ Single velocity minimal 
STOL navigation. navigation error. 
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Operations rate per hr 


Minimal spaaing— Finally, consider the case (fig. 9) in which the spacing 
is reduced to correspond to the increased navigation accuracy. To determine 
the extreme situation, a minimal spacing value of 0.1 km was used. With this 
spacing, the landing rate is determined almost exclusively by the runway occu- 
pancy time and, to a lesser degree, by the error in touchdown velocity. The 
addition of STOL traffic increases the landing rate, and a STOL-only runway 
again has a higher landing capacity than the CTOL-only runway. Moving the 
STOL merge point closer also increases the landing rate by decreasing the 
allowance required to compensate for the uncertainty in final velocity. For 
this case, times are available for only a few takeoffs, and then only with 
mixed traffic, when the distant merge points are used. 


Dual-Velocity Approach 

For a dual-velocity approach, all CTOL aircraft use the same constant 
speed from the CTOL merge point until time to change velocity, about 5 or 6 km 
from the threshold. As a result, the location of the merge point for CTOL 
vehicles has very little effect on the landing or takeoff rate; therefore, it 
is not shown as a parameter in the figures. 

The landing rates and the total operation rates for the dual-velocity 
approach are shown in figures 10 through 15. 



Figure 9.— Single velocity 
minimal spacing. 


Figure 10.— Dual velocity 
present standards. 
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Present standards— The results for this case are shown in figure. 10, and 
the trends are generally the same as for the constant velocity approach. 

Adding STOL to the traffic mix reduces the landing rate but increases the 
total operation rate. The landing rate is controlled strictly by airborne 
spacing, usually at touchdown but, in some cases, at the velocity change 
point. 

Nominal— The reduced spacing for STOL vehicles (fig. 11) now results in 
a strictly increasing landing rate as the proportion of STOL vehicles 
increases and in an increasing total operation rate. The controlling factor 
is the airborne spacing for a CTOL first aircraft and either airborne spacing 
or runway clearance, depending on WspoL* f° r a STOL first aircraft. The land- 
ing rate for this case is somewhat larger than for the corresponding single- 
velocity case (fig. 5), but the total operation rate is less. 

Reduced runway time— When the runway occupancy time was reduced (fig. 12), 
the landing rate increased compared to the nominal (fig. 11), for Wg T Q^ = 7 km. 



Figure 11.— Dual velocity Figure 12.— Dual velocity reduced 

nominal case. runway time. 
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where runway occupancy was the controlling factor, but otherwise was unaf- 
fected. The total operation rate, however, increased considerably compared to 

the nominal. For this case, increasing the STOL traffic increased the landing 

rate and also increased the total operation rate. 

Improved STOL navigation— Again, as in the single-velocity case, improv- 
ing the STOL navigation and guidance (fig. 13) had the same effect on the 

landing rate as reducing the runway occupancy time (fig. 12). The total 

operation rate, however, is essentially the same as for. the nominal case 
(fig. 11). Both landing rate and total operation rate increase with increas- 
ing STOL traffic. 

Minimal navigation errors— Further reductions in the navigation and 
guidance errors (fig. 14) had no effect on either the landing rate or the 
total operation rate. 



Figure 13. Dual velocity improved Figure 1 14. ~ Dual velocity minimal 

STOL navigation. j navigation errors. 
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Minimal spacing— Reducing the 
spacing to a minimal value (fig. 15), 
however, increased the landing rate 
dramatically over the nominal, but 
left no time for takeoffs. 

Comparison with single-velocity 
case— Generally, the landing rate for 
the dual-velocity case is somewhat 
higher than that for the correspond- 
ing single-velocity case. This 
higher landing rate tends to accom- 
pany a somewhat lower total opera- 
tion rate. 

The trends for both the dual- 
velocity and single-velocity cases 
are very similar with changes in 
spacing, navigation, and runway occu- 
pancy time. However, as the percen- 
tage of STOL traffic increases, the 
landing rate always increases in the 
dual-velocity case; in the single- 
velocity case, the landing rate first 
dips and then increases. 


CONCLUSIONS 


The following conclusions are 
based on the philosophy that landings have priority over takeoffs. 

Adding STOL traffic to a CTOL runway causes only a minimal decrease in 
the landing rate for a constant-velocity approach, while for the dual-velocity 
approach, there is an actual increase in the landing rate. This assumes that 
the STOL navigation is sufficiently better than existing CTOL navigation so 
that the STOL spacing can be reduced from 6 to 2 km (3 to 1 mile)-. With this 
spacing, the landing rate on a STOL- only runway will be greater than on exist- 
ing CTOL runways . 

Under all circumstances, including the use of present standards, the 
addition of STOL traffic to a CTOL runway will allow an increase in the total 
operation rate (i.e., landings plus takeoffs). Thus, with mixed STOL-CTOL 
operations, there is frequently time for a takeoff without disturbing the 
landing sequence so that mixing takeoffs and landings on the same runway is 
feasible and perhaps desirable. 

To increase landing rates (from present values) , airborne spacing must be 
reduced, which, in turn, requires improvements in navigational accuracy. 



50 

STOL, percent 


Figure 15.— Dual velocity 
minimal spacing. 
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Reducing runway occupancy time has little effect on the landing rate but does 
allow for an increased total operation rate. 

Generally, the landing fate is higher if the STOL merge point is as close 
as possible to the runway threshold, thus taking maximum advantage of the STOL 
maneuvering capability. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, October 29, 1973 


- -w — 

i 

:3vvTr 







21 



APPENDIX A 


DERIVATION OF EQUATIONS FOR T l2 


The kinematic equation of motion of an aircraft approaching the runway 
is given by equation (2), assuming that the aircraft touches down at t = 0. 
This is true for the first aircraft, but the second aircraft touches down at 
time T \2 after the first.. The equations of motion for the two aircraft are 
then: 


*1 




(Al) 


/ \ 



-%(* - 

*2 ^ Da z 

x 2 = • 

CM 




x 2 > V 
2 


(A2) 


From these equations, one can calculate the times 
at specific points. These times, and their equations, 
first aircraft is at its merge point, mi. 




m i 


when the aircraft are 
are: the time the ! 



i 

I 



(A3) 


the time the first aircraft is at the merge point for the second aircraft, 

m 2 . 



-m 2 /V, 


m 2 


D. 


a , 


v V. 


a , 


V. 


a , 



m 2< D a \ 


m 2 > D ai ) 


(A4) 
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and, finally, the time at which the second aircraft changes its speed. 


Z? 


a , 


T a _ = T 


12 V. 


(A5) 


b 2 


The equation that defines the reserved airspace for the first aircraft is 
obtained by combining equations (4) and (5) : 


Si («i) Jfi 


Xi = 


( 


+ 1 + 


S °1 ' S fl 


2m i 




fl 


Ki%i , Xi < mi 


(A6) 


where 


- S 


K : = 


( 

\ + 2 mj 


A' 


Note that since no tests are needed at times before the first aircraft arrives 
at its merge point, a:]! is not needed for > mi. Similarly, the reserved 
airspace for the second aircraft is 


S 2 (# 2 ) 


x 2 = x 2' 


= < 


% 


S 0 2 " ^ 


+ 1 - 


X 2 


2m n 
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x 2 < m 2 


x 2 > m 2 


Sf 

- + ZZ 2 x 2 , xi < m 2 


+ x 2 , x 2 > m 2 


where 


*2 = 
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Assume that S Q >_ Sf, that is, the required separation at touchdown may be 
less than (or equal to) that at the outer marker, but never greater; thus, 

K\ > 0. However, if 



> m 2 


! > 



then K 2 < 0. Consider first the situation when K 2 > 0 (fig. 16). When 
'x 2 :> ~m 2) the separation is S Q / 2; when x 2 < m 2 , the separation is 
(S--/ 2) + K 2 x 2 . Next, figure 17 shows the situation for K 2 < 0, when the 
12 


Flight path 


m 2 



Figure 16.— Reserved space before 
second aircraft, K 2 > 0. 



Figure 17.—' Reserved space before 
second aircraft, K 2 < 0. 


first aircraft will land before the second aircraft reaches its merge point, 
so that airborne separation after the merge point is not required. Thus, for 
K 2 < 0, the second part of the abqve equation for x 2 should be used. With 
these restrictions. 



fz 

- + K 2 X 2 , 

x 2 < m 2 

and K 2 > 0 j 

x 2 ' = • 

s 

O 2 

■ T~ + X 2 

X 2 > ™2 

or K 2 < 0 I 


(A7) 


With the equation of the boundaries, the next step is to calculate T X1 
so that the boundaries do not cross but touch at some point. This will 
clearly determine the minimum T\ 2 since, if it were smaller, the separation 
criteria would not be met at some points along the trajectory. Thus 
(t) - x 2 {t ) = 0, is solved at selected values of time, for T l2 . This 
equation must be solved for several different values of t, namely, at the 
final point, t - 0, at the merge point, t = max(Tm a >T mi ) , and for the case 
when the second aircraft slows down during the approach, at the velocity 
change point t = T a% . To minimize confusion, the time difference from each 
equation is labeled x, and !Fi2 max ( T f , ^n , T med> r o) - (The value o t g 
is determined later.) Thus, the following three cases are obtaine . 
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I: 

IIA 

IIB 

III 


x](0) - aj 2 '(0)"= 0 , solve for 

*1 Wmc? ' ^'(^J = 0 , ^2 < ‘ 

x{(T ml ) - x 2 ' (T,,) = 0 , m 2 > 


ma 


solve for t 


tn 


*i(s) " ^'(S) = 0 • m{ - T ma tT mO K T a z K 0 ’ solve for \ ed ) 


(A8) 


If the restriction for case III test is not met, then a test for case III 
is not needed; for convenience, we set T mg j = 0 so that it cannot be' selected 
as the maximum in determining T \ 2 . To solve equation (A8) for x, use equa- 
tions (Al) through (A7) . Since these are conditional equations, depending on 
the relationships between x 2 {t) , D a 2 , m 2 , etc., many subcases arise. Each 
must be solved for and then checked to ensure that all conditional relation- 
ships are met, and, finally, the defining conditions for each subcase must be 
determined in terms of the input quantities, D a , Sj>^ , V a ^ , F^, and 

i = 1,2. This procedure is followed for each subcase in the following sec- 
tions. A small sketch shows position versus time and indicates the location 
of the input quantities for each subcase. 


CASE I 


In this case, t = 0, so the following equations are needed 
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+ K 2 x 2 , x 2 < m 2 and K 2 > 0 


1 


~Y~ + x 2 , x 2 > m 2 or K 2 < 0 
These equations lead to four subcases, each of which is considered separately. 
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(A10) 


(All) 
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Case la: x 2 1 D a and (x 2 < m 2 anc ^ ^2 - 0) 



Combining the appropriate 
equations from (A9) to (All) gives 


"f ! + ^2 


/in 2£ 2 7^ 


If this value for ri 2 is used 
in equation (A10), 


Sketch (c) 


S fl + ^2 

* 2 ~ 2 Ko 


To meet the conditions for this case. 


x 2 < D 0 


x z < m 2 


S„ + S„ 

f \ ?2 

2K 2 

S f + S f 

J 1 J 2 

2 K 0 




m 2 


K 2 > 0 

These three tests can be combined into two, since K 2 > 0, if both sides are 
multiplied by K 2 : 

i S f 1 + S f 

x 2 5 D a'\ 2 ~ K2 - D a, 


I S f +S f 

I J 1 J 2 

11 ■ 1 x 2 < m 2 | *■ i 2 1 -^2^2 

which are automatically not satisfied if K 2 < 0. 
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Case lb: x 2 > D a and (x 2 1 m 2 and K z > 0) (sketch (d)) 

If these conditions are used in 
equations (A.9) to (All) , 

V n 


T = 


% * S f 2 ° a 2 / '°2 

1 - 


fin 2 K 9 V. 


2 a„ 
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12 ‘ Si 1 " Vi I 
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a . , 
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Sketch (d) 


+ V 


a , 


' s fl * s f 2 \\( 1 _ S 


2K 2 V 


a r 


V 


<-2 a 2' V u 2i 

To meet the necessary conditions, with K 2 2- 0, one must have 
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S h * S /2 

2K 2 


x 2 > D { 
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2 > 


S-p + S p 

f 1 ?2 


x 2 < m 2 


K 2 m 2 


Case Ic: x 2 < D n and ( x 2 > m 2 or K 2 < 0) (sketch (e)) 

u 2 

Again, solving equations (A9) to 
(All) gives 


S f l * S 0 Z 
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fin 2V b 


S f , * V 
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Sketch (e) 
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The condition tests are 



Sketch (f) / "J a z\ S fl + S °2 

" Da 2\ \) = — ~ 


The condition tests are 






x 2 > rn 2 




Remarks ' 

Although all four cases have been discussed, the tests do not appear to 
guarantee that all possible conditions were included. However, if the value 
of K 2 is substituted into the second test for case la, one obtains 
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or 


S f ! + S f 2 


2 




< m 2 


which is the identical test used in 
cases Ic and Id, thus showing that 
all possible conditions were included 
Figure 18 is a flow chart of the case 
I tests. The test used in the figure 
is that described under cases la and 
lb. 


. CASE IIA 



Here, the minijnum separation is 
determined when the' 1 first aircraft is 
at the merge point. The merge point 
is assumed to be m 2 , so that Xi = m 2 , m 2 < and t 
needed obtained from equations (Al) to (A7) , are 


Figure 18.— Calculation of x 
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Note that only one equation for x 2 is needed here since x 2 > aa^ = m 2 : Further, 

since K 2 does not enter the equation, its value need not be considered. 
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For each of the four cases here, solve equations (A12) to (A15) for T i 2 ~ x i n • 
The solutions and the required tests are as follows: 


Case IIAa: X\ < D a , 



Sketch (g) 


x 2 - D a 2 ( sketch Cg) ) 
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Case IIAb: 
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Sketch (h) 
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Note: D o2 shown <m2, but moy be >n>2 

*1 > D a 1 m 2 > D a x Sketch (j) 

S ?l + S ° 2 

x 2 > D a 2 2 + KlTnz > °a 2 

Figure 19 is a flow chart of the case IIA tests, A4. 



Figure 19. 
CASE I IB 


Again, the minimum separation is determined when the first aircraft is at 
the merge point. Here, however, the merge point is at so x\ = m± and 
m 2 > mi, with t = T m . The equations needed are 
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Here three conditions must be considered, which lead to eight cases. For each 
case solve equations A(16) to A(19) for T 12 = t . . 


Case IIBa:' x 2 < D a , (x 2 < m 2 and K 2 > 0) , x 1 < D a ^ (sketch (k)) 



S-p + S-p- 
J 1 •’2 m i 




T in ’ Z 2 F &2 
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Sj? + Sy 

x 2 1 D ~ — + K l m l tWa, 


Sketch (k) 

Again, as in case I, since K 2 - 0, both sides have been multiplied by K 2 . 
The resulting test will fail if K 2 <0: 
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X 2 - m 2 


+ Kim i < K 2 m 2 
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Case IIBc: .. x 2 < D a ^, (x 2 > ^2 or ^2 < °) > X 1 < ^ (sketch (m)) 



x 2 > m 2 -y —2 + *l m l > «2 Sketch (m) 


Note that, as in case I, substituting the value of K 2 into the test for 
x 2 < m 2 of case IIBa or IIBb gives the tests determined above. 

x i ^ \ 1 D ai 
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Case IIBd: x 2 > D Q , {x z > m 2 or K 2 < 0) , x 1 < D ( 
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Case IIBh: 


x r > D 


a , 


(a? 2 > m 2 or ^2 < 


0 ), *! > D a 


1 


(sketch (r)) 



Figure 20 is a flow chart of case I IB. 


CASE III 


For this case, the minimum separation is determined when the second 
aircraft is at t = T a 2 . Note, however, that this test applies only if 
Max (T m a , ^m\) < T a2 < 0. Since the value of T a depends on the value cal- 
culated from T\i - t med> the procedure is to first calculate % ec f, according 
to the conditions imposed by x\, x-i , etc., and then, for these conditions, 
check that the value of T a 2 is acceptable. The equations used are again (Al) 
to (A7) , with t = T a : 


= < 


~V h T„ 

b l a 2 


V , 


a , 


, D a 1 1 

V a l 


Vv 


X 1 £ \ 


W > X 1 > D a 


l / 


(A20) 



(A21) 
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Figure 20.— Calculation of x . for m 7 < m, . 
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x 2 = 


+ K 2 x 2 , x 2 < m 2 and K 2 > 0 


+ x 2 , x 2 > m 2 or K 2 < 0 


(A23) 


D 


a , 


r a 0 = ^12 


Vi, 


GA24) 


which leads to four subcases. Further, to test the acceptability of the 
value of T a ^> one must also consider the equations for T mi if m 2 > m i , and 
T a if m 2 ± mi . These equations are 
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T = <. 
ma 


V b , 




m 2 a 1 


XV V 

' a l a l 


/ m. 


T = < 
m l ' 




m 


D, 


a. 


7. 


a. 


V-L 


CL, 
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' v “, r «, 




m 2 < 


> D 


a 


1/ 


m 2 < mi 


m l < 0, 


< 2 , 


> D a, 


m 2 > mj 


(A25) 


Note that the form of these equations is identical, the only change being in 
the choice of mi or m 2 . In the sequel, these equations will be combined into 
one equation, with m . = minG^i, m z~) • 
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Case Ilia: (x 2 <m z and K 2 >.0), x± <D ai (sketch (s)) 
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a? 


Sketch (s) 


To check for applicability, use 


+ S fz 


T a 2 " 0 


< K l D a, 


Note that this will not be satisfied unless K 2 > 0, so a separate check for 
this is unnecessary: 
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, rn . > D a 
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Note that when %iin > i> the left side of the last part of the equation is 
positive by the test for xi <D ai , and the right side is negative since 
m min > D ai’ Therefore, the inequality is automatically satisfied if m m i n > ^a 1 
and need not be checked. 
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Sf + ,%* : 
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Note that K 2 does not enter these equations so its sign is not important: 
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where the second test is automatically satisfied as in case Ilia. 


Case IIIc: 


(r 2 < rn 2 and K 2 >0), 
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This test is satisfied since the test for x^ > D a is satisfied, and therefore 
need not be considered: 1 


(*) 


max 


( Tm a’ Tm l) < \ 


S jy + S jy 

J 1 1 2 


Kim . > K 2 D 

■L mi -n *- l 


min 


a. 


m . < 

min 




1 


m . > 

mm 



Again, as in the previous two cases, the first of these two tests is included 
in the previous tests, but now the test can never be satisfied; therefore, it 
need not be tested. Also, if the second test is satisfied, then m m x n > ^a\> 
if this test is not satisfied, case III is not needed (so tests for m . > D a 
are also not needed) . 1 
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This test will always be satisfied since xi > T> a and therefore need not be 
checked: 1 
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As in case IIIc, the first part of this test can never be satisfied and so 
need not be tested. Also, the test for m m ^ n > D al can be omitted. As shown 
in figure 21, where the tests for T a < 0 or maxCr^a* < T a 2 are not sat " 

isfied, the value of is set to zero, so that it is ignored in the maximi- 

zation process for T ^ 2 . 



Figure 21 .— Calculation of t mg 


If T\2 = max (xi n , tfin> T meS > adequate separation will be maintained 
while the aircraft are airborne, but the second aircraft may still land before 
the first aircraft is off the runway. To prevent this, an .additional quantity, 
t 0 , is computed which controls the spacing between the aircraft so that the 
first aircraft will be off the runway before the second lands despite errors 
in the aircraft's maneuvers. 


42 


















RUNWAY SEPARATION 


The interarrival time required for runway separation is defined by 
equations (10) and (11) . The partial derivatives used in equation (10) are 
defined as 


3Tl2 

got 

30 . 


3 T 

nom 

30 . 
x . 


where T is the negative of equation (3) or (A3) : 

( t 



m < D a 


m > D, 


a 


With this equation, along with equation (1), the various derivatives can be 
calculated as follows: , 
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All the quantities needed to calculate in equation (12) are now 

available. rU)m 
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PROBABILITY OF ABORT 


In equation (10) , a value for A T is determined. This quantity is an 
allowance to be added to the runway occupancy time to reduce the probability 
that errors in arrival position and velocity will cause the second aircraft 
to abort, that is, if the first aircraft is not off the runway when the second 
aircraft wants to touch down. Thus, it is desirable to determine analytically 
the probability of abort so that a desirable value for K is obtained. 

The probability of abort is a function of four random variables (6V a , 

&Vb> dp, and 6T v f~) for each aircraft (a total of eight variables), each with 
an assumed triangular distribution. The resultant distribution of deviations 
of interarrival times (sketch (w)) is essentially Gaussian. Also shown is 



AT 


Sketch (w) 

AT, K times the. root sum square of half the ranges of the n individual tri- 
angular error sources. The probability of an abort is shown by the shaded 
area. To estimate the probability of abort, first assume that all variables 
have been normalized and have uniform ranges. Also, note that a triangular 
distribution function is the convolution of two identical uniform distributions. 
Therefore, the probability of abort of the n triangular distributions is the 
same as that of 2 n uniform distributions. Thus, given the density function 
for 2 n uniform identical random variables. 


4 00 = { 1 ’ 
X l l 0 , 


o < a < 1 
else 


(A26) 


determine the probability that 


2 n 



i= l 


where T s = ^n. Thus, T s corresponds to AT, with K = 1 . Using Laplace trans- 
forms, the moment generating function <J> for a single uniform variable is 

v 
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From, this, the overall function <b is 

3 



(A28) 


Taking the inverse Laplace transform gives the density function. 
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which can be integrated to give the probability distribution function: 
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CA30) 


To determine the probability of abort, evaluate equation (A30) at t = x - T s 
where x = n is defined as the mean of 2 n identical uniform distributions, that 
is, the probability that Z < x - T g is F z (n - Jn ) . This evaluation yields 
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F {n - /ri) = 


(-D 

k\{2n - k ) ! 


, k > n - Jn 


j— r - 

[n - vn - k) , k < n - vn 


max 

■E* 


f-i) /- 2n 

(2 n - k ) ! (n ' /n " 


where fc m ax ; *- s t ^ e largest integer k so that k < n - Jn or k max = int[n - /n] . 
Table 3 gives, for various values of n, the values of ^ , F z (n- /n) , 


TABLE 3.- PROBABILITY OF ABORT FOR VARIOUS n 


n 

k 

max 

F z [n - •In') 


1 

0 

0 

“oo 

2 

0 

.00491 

-2.58 

4 

1 

.00615 

-2.50 

8 

5 

.00670 

-2.47 


and Zg, the normalized Gaussian variable with the same area under the left- 
hand tail as F z {n - Jyl) . 


Figure 22 is a plot of Zg versus 
n. The previous calculations assume 
that all error sources are identical, 
that is, each has an identical effect 
on the landing time whereas, in fact, 
each has a different effect. To 
consider nonidentical sources, assume 
that n - 1 sources had identical 
effects and that the effect of the 
nth source was to vary from that of 
the others to zero. As it varies, 
the width of the probability curve 
decreases, and, by definition, the 
value of A T also decreases. The net 
effect should be that the area to the 
left of AT will change monotonically 
from the probability of abort due to 
n sources, to that due to n - 1 
sources. From figure 22, note that 
this probability is very insensitive 



n 


Figure 22.— Effect of number of error 
sources in probability of abort. 
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to n and, furthermore, the probability of abort decreases as n decreases. 
Therefore, the probability that the actual error will exceed the RSS of 1/2 
the ranges (which implies K = 1) occurs with a probability of less than 0.006, 
as stated in the text. The Monte Carlo runs showed the probability of aborts 
from 0 to 0.006, which agrees with the above analysis. 
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APPENDIX B • 

PROGRAM LISTING AND DESCRIPTION 


A listing of the complete program, consisting of a main program and two 
subroutines,- HEAD 1 and OUT, is presented. The main program first calls HEAD1, 
which reads input data, sets up the calculations, and prints heading informa- 
tion. The main program then does the major portions of the calculations and 
some output, and finally calls OUT for the remainder of the output. The pro- 
grams have many comment cards and it is hoped that they are self-explanatory. 

The computations for a single case are done from a data array called 
ACDTA, and the purpose of HEAD1 is to change the data in ACDTA as required. 
This allows a kind of DO LOOP action, under the control of HEAD1, so that a 
variety of conditions can be run simply. A listing of a sample input deck is 
discussed in detail to show how the required, inputs are made. 


INPUT DECK 


1 VEL VARIOUS i 


Mixen stol ctul euuAi spacing 1 


•35' *30. . ' 40 1 ,19., 
•70' 


The inputs are free form as controlled by the Ames INPUT routine. 1 This 
routine continues to read data cards and put the data into the named locations 

until a * is read. At that point, 
the program proceeds to the next oper 
ation, and subsequent data cards are 
read by a subsequent call to INPUT. 

Figure 23 shows the sample input 
deck, where card 1 is a title -card 
and card 2 defines the number of land 
ings. to be used in the statistical 
portion of the program. Cards 3 
through 8 give information concerning 
the aircraft, one card for each air- 
craft type. Up to 12 aircraft types 
can be defined and loaded into AC1 
through AC 12, which are rows of the 
ACDTA table. The variables for each 


AC 1=1 9. ,6. ,6. ,60., 60., 60. ,90. ,2. ,2.5,2.5,30. , 
AC?=19.,6. ,6. ,60., 65. ,65., 90., 2. ,2.5,2.5,30., 
AC3=19. ,6. ,6. ,60. ,70. ,7 0. ,90. ,2. ,2.5 ,2.5,30. , 
AC 7= 2. , 2 . ,2 . » 30 . , 30. , 30. , 60. , 2 . , 2 . 5 . 2 . 5 , 1 0 . , • 
AC8 = 2.,2.,2. ,30. ,35. ,35 . , 60. , 2. , 2.5 , 2 . 5 , 10. , • 
AC9= 2. ,2. ,2. ,30. ,40. ,40. ,60. ,2. ,2. 5, 2. 5, 10., 
'ACPROBs • 30 
C HNGAC = • SO 1 
TVF=900.,. 

SF-.l, 

NORUN= 1 ’ , - . 

CHNGAC* * 30* , *35', '40' , 

TVF=9O0. , : 

SF«.l, 

M=2. ,7 • »• 

ACPROfi= • 30* ,120.', '35* ,90. , '40* ,57. , 

CHNGAC =• 60' , • 65' , '70* , 

M= 7 . » 

N0RUN= 1 , 


1 65' ,30.. ' 70', 19. 


CHNGAC = * 30’ 
M= 2 . » 7 . , 
ACPROR= * 30* 

CHNGAC- • 60' 
M=19., 
N0R1IN= 1 , 


85 

0- 


13 

15 

1 6 
17 
13 
19 


24 

25 


27 

23 


•35 ’ ,30. , '40* , 19. , 


40. , ' 65* , 30., *70' ,19. 


CHNGAC- 
M=2. ,7.. 


3 s * 

35 

36 

37 


V, 


a> 


Figure 23.— Input deck. 

various aircraft types to be considered, 
new proportion card is entered. In this 
are to be used. 


aircraft type are m, S Q , Sf, T c 
Vb> T vf> R&P> R&V a> R&Vb, R&T v f, and 
type (in units of km, sec, and m/sec) 
Card 9 defines the proportion of the 
and this proportion is used until a 
case, only the three STOL vehicles 


^his allows free-form input similar to NAMELIST, and a description of 
the subroutine can be obtained from the Ames Computer Library. . 
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To reduce the number of input cards required to run multiple cases, the 
input procedure was arranged to set up a series of DO LOOPS. With this pro- 
cedure, it is necessary to first indicate which of the aircraft types is to 
have variables changed by the DO LOOPS and then to indicate what variables are 
to be changed, and their new values. If more than one value is specified for 
a variable, the program will compute for all values specified, as in a DO LOOP 
If more than one variable is to be changed, nested DO LOOPS are created. Thes 
DO LOOPS remain in force until a new indication of aircraft types is specified 

The particular set of input data shown in figure 23 computes landing 
rates for the following cases: 


STOL only 


W ST0L " 2 5 

7 



75 

percent 

STOL 

W ST0L = 2 ' 

7 

W CT0L = 

19 

75 

percent 

STOL 

m ST0L = 2 ’ 

7 

W CT0L = 

7 

50 

percent 

STOL 

m ST0L = 2> 

7 

m cm = 

7 

50 

percent 

STOL 

m ST0L = 2j 

7 

m CT0L = 

19 


The remainder of the input cards provides the information to run these cases. 

Card 10 indicates that, for the first case, some data for the three 
aircraft types indicated (CTOL) are to be changed and cards 11 and 12 indicate 
the required changes to T v f and Sf. Card 13 indicates that no calculations 
are to be made at this point (since more changes will come in the next call to 
INPUT). Cards 15, 16, and 17 make the same changes for the STOL vehicles, 
while card 18 indicates that computations should be done for rngjoL = 2 and 
m ST0L = 7 and has thus set up a one- level DO LOOP that will be gone through 
twice. Up to four values of m could have been provided, and additional values 
for T v f and Sf could also have been provided. If two values had been speci- 
fied for both T v f and Sf, then there would have been a three level DO LOOP, 
with two values for each level, resulting in eight sets of computations. 
Quantities that can be changed in this fashion and used to control DO LOOPS 
are m, S Q , Sf, T 0 , V a , V^, T v f, E6p, R6V a -, R&Vb, and R6T v f. If all these vari 
ables are used, an eleven-level DO LOOP will exist, with up to four values 
for each level, resulting in 10 7 cases. Obviously, care must be exercised 
here. 


Card 19 again indicates the end of a data read, and only two cases (STOL 
only, msTOL = 2, 7) will be computed. The output from these two cases is 
shown in figures 24(a) and 24(b). Card 20 presents a new set of proportions 
of aircraft types and, since there is no reference to CHNGAC, the previously 
set-up DO LOOP will be used for the next cases (75 percent STOL, wgjoL = 2, 7, 
w CT0L = 19) • The output from the first of these cases is shown in 
figure 24(c). 

Cards 22 and 23 indicate a change to the CTOL aircraft data, to set 
m CT0L = 7, and cancel the initial DO LOOP setup. Card 24 indicates no compu- 
tations. Cards 26 and 27 change variables for the STOL aircraft, and again 
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set up the one-level, two-value, DO LOOP (75 percent STOL, msxOL = 2, 7, 
m CTOL = 2) . Card 29 specifies 50 percent STOL, running the same DO LOOP 
(50 percent STOL, msTOL = m CT0L = 7) . On card 37, the last case has 

been set up to run (50 percent STOL, msTOL = .2, 7, WcTOL = 19 ) • The program 
will continue to read data cards until they are gone, and then it will stop. 
This particular input deck will compute data for 10 cases (as previously 
stated) . ' 

Two additional subroutines are called by the program. The first, 
RANDU(IX,NX,R) , is a random number generator and (as used here) will return 
a random number, R, with a uniform probability distribution from 0 to 1 . The 
variable IX is a starting point and NX is the starting point for the next 
random number. The second routine, VTITLE, prints a title line and page 
number at the top of each page of output. The title is read from the first 
data card. 


DESCRIPTION OF OUTPUT 


The first three pages of output from the sample input deck are shown in 
figures 24(a) through (c) . The output includes not only the information dis- 
cussed in the text, but also some material which is not discussed. This 
description defines all of the data in the output. 

In figures 24 (a) -(c), the first line is the case title and the second two 
lines are general heading information. This is followed by the "A/C DATA 
TABLE," which gives all the input quantities for each type of aircraft used 
for this particular case. In addition, the individual partial derivatives 
required in equation (10), along with the root sum square, are given in the 
last five columns. 

The next two lines define the probability that the next aircraft is of a 
specified type. This is followed by a table defining the interarrival times 
and available takeoff times for all possible sequences of aircraft. The first 
line defines the type of the second aircraft, and the first column defines the 
type of the first aircraft. For each possible sequence, four entries are made 
in the table*. The first, PROB, is the probability of that sequence, obtained 
by multiplying the probability of the individual types. The second item, T12, 
is the required interarrival time, as calculated by equation (12) . Immedi- 
ately following this is a letter, which defines the x with the largest 
value, using the code indicated. The last item is the time available for 
takeoffs, which is the interarrival time minus the runway occupancy time of 
the first aircraft. 

Some histogram information is presented next. The total number of 
aircraftin the histogram is specified. The histogram shows the distribution 
of interarrival time and available takeoff times, as calculated from the 
information immediately above. Also, average, maximum, and minimum times are 
shown. 
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Figure 24(a) - Sample output. 
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Figure 24(b).- Sample output - continued 
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Figure 24(c) - Sample output - concluded 



Next, data on the landing rate are presented. The Monte Carlo process 
was used to determine how many of the landings were required to abort because 
an aircraft was attempting to touch down before the runway was clear. The 
abort ratio (number of aborts/total number of landings) is printed, along with 
the nominal (or planned) number of landings and landing rate. The actual land- 
ings referred to in the last two columns are the nominal landings minus the 
aborts . 

Finally, takeoff information is presented. The first line of this table 
(fig. 24(a)) is read as follows: if all aircraft use 36 sec for takeoff, then 

there will be sufficient time for 1089 aircraft to take off without interfer- 
ing with the landing vehicles. Since there are 302 landings being considered 
(given on the previous line), there will be a total of 1391 operations, at a 
rate of 90.7 operations/hr. The last two columns consider the effect of 
aborts on the number of operations and the operation rate. Since there are 
no aborts (see previous line), the actual and nominal operations are the same. 
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